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We investigate the different regimes that emerge when a periodic driving force, the rocking force, 
acts on a collective feedback flashing ratchet. The interplay of the rocking and the feedback control 
gives a rich dynamics with different regimes presenting several unexpected features. In particular, we 
show that for both the one-particle ratchet and the collective version of the ratchet an appropriate 
rocking increases the flux. This mechanism gives the maximum flux that has been achieved in a 
ratchet device without an a priori bias. 
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I. INTRODUCTION 

Ratchets can be viewed as controllers that act on 
stochastic systems with the aim of inducing directed mo- 
tion by breaking of thermal equilibrium and certain time- 
space symmetries [l[ . As usual in control theory 0] , these 
systems are divided into open-loop ratchets when the 
actuation does not use any knowledge of the state of the 
system; and closed-loop ratchets 3, 4], when information 
on the state of the system is used to decide how to op- 
erate on the system. These closed-loop ratchets — also 
called feedback or information ratchets — have recently 
attracted attention as Maxwell's demon devices that are 
capable of maximizing the performance of ratchets ||. 
They may also be relevant to get insight into the motion 
of linear, two headed, processive molecular motors Q. In 
addition, experimental realizations of feedback ratchets 
have been recently proposed [1, 0] and implemented @ 
due to their potential relevance as nanotechnological de- 
vices. 

A relevant class of ratchets are flashing ratchets, which 
operate switching on and off a spatially periodic poten- 
tial. Flashing ratchets have been studied in both open- 
loop (e.g. [Eil) and closed-loop (e.g. 0,13, E3]) schemes. 
A generalization of these ratchets are pulsated ratch- 
ets [l|, in which the amplitude of the ratchet potential 
is modulated in time, but not necessarily flashed on and 
off. On the other hand, rocking ratchets operate thanks 



to a periodic driving force, and thus they perform an 
open-loop control. Rocking ratchets reveals a rich dy- 
namics, which includes current reversals, distinct stable 
trajectories, and quantization of the deterministic cur- 
rent [HI, HH, The combination of open- loop pul- 
sated ratchets and rocking ratchets has been studied in 
Refs. (l4j |. giving the possibility of a reverse of the sign 
of the flux with respect to the simple rocked ratchet. 
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In the present paper we study the effects of adding 
a periodic driving force that rocks a feedback-controlled 
flashing ratchet. We analyze the new intriguing dynam- 
ics that emerge due to the interplay between the feedback 
control and the rocking. In particular, we show that the 
rocking of a feedback ratchet allows the system to im- 
prove the flux performance. The optimization of the flux 
performance of ratchets is potentially relevant for their 
nanotechnological applications, and the enhancement of 
the flux performance in flashing ratchets due to feedback 
has been recently verified experimentally Q. We show 
here how this flux performance can be further improved 
thanks to the effects produced by an additional rock- 
ing force. In the next section we describe the rocked 
feedback-controlled ratchet, and after, in Sec. IIII1 we 
study the one-particle ratchet. The collective version of 
the ratchet compounded of more than one particle is ana- 
lyzed in Sec. lIVl in the regimes of few and many particles. 
We finally summarize and comment our main results in 
Sec.lV] 
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II. ROCKED FEEDBACK-CONTROLLED 
RATCHET 

Let us consider N Brownian particles at temperature 
T in a periodic potential V{x), the ratchet potential. The 
state of the system is described by the positions Xi(t) of 
the particles (i = l,...,iV) satisfying the overdamped 
Langevin equations with a fluctuating (rocking) force of 
amplitude A and frequency il, 

j±i(t) = a(xi(t), . . . , x N (t), t)F(xi(t))+AcoB(nt)+ti(t). 

(1) 

Here, F(x) = — V'(x), 7 is the friction coefficient (related 
to the diffusion coefficient D through Einstein's relation 
D = ksT/^f), and £i(t) are Gaussian white noises of zero 
mean and variance (£,i(t)^j(t')) — 2jkBTSijS(t—t'). Note 
that the control parameter a depends explicitly on the 
state of the system. Therefore, this ratchet is feedback 
controlled, what implies an effective coupling between the 
particles. 

In order to quantify the inducted directed motion a 
relevant quantity is the stationary center-of-mass velocity 
or flux defined as 



lim 
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( Xi {t) - Xl (0)) 
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Due to the long-time limit and the average over real- 
izations this asymptotic center-of-mass velocity does not 
depend neither on the phase of the fluctuating force nor 
on the initial particle positions [l[ . We shall consider the 
relevant control policy that maximizes the instant center- 
of-mass velocity introduced in [1]. In this feedback pro- 
tocol, the controller computes the force per particle due 
to the ratchet potential if it were on, 



1 N 

f( Xl (t),...,x N (t)) = -Y. F ( x ^ 



(3) 



and switches the potential on (a = 1) if f(t) is positive or 
switches the potential off (a = 0) otherwise. Therefore, 
the feedback control protocol considered is 

a(x 1 {t),...,x N (t))=Q(f(x 1 (t),...,x N {t))), (4) 

with 9 the Heaviside function [Q(x) = 1 if x > 0, else 

e(x) = o]. 

The graphs that illustrate the results of this paper have 
been obtained considering the periodic asymmetric po- 
tential 



V(x) = 



2V 



sin 



2ttx \ 1 . / 4ttx \ 

t) + 2-(t) 



(5) 



which has potential height Vq and period L; see Fig. [T] 
We can introduce an asymmetry parameter a for the po- 
tential such that aL is defined as the distance between a 
minimum of the potential and the first maximum at the 
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FIG. 1: (Color online) Ratchet potential V (solid red line) 
[Eq. ([5} with Vo = 5], one-particle effective potential V e s 
(dashed green line), and one-particle periodic effective po- 
tential V c g Cr (dotted blue line). Units: fcsT = 1, L = 1. 



righthand-sidc. The potential in Eq. ([5]) has an asym- 
metry parameter of a = 1/3. We have also performed 
computations with other potentials and found analogous 
results. 

In order to numerically compute the flux (J2J) we have 
performed numerical simulations of the Langevin equa- 
tions (p} [with the control parameter a given by Eq. f4|] 
by using the Euler-Maruyama algorithm [l5j| . We have 
verified that our numerical results were not affected by 
systematic errors due to time discretization, initial tran- 
sitories or finite number of realizations. The characteris- 
tic time that takes the system to diffuse the distance aL 
of the uphill part of the ratchet potential is a 2 L 2 / (2D) , 
while the characteristic time that takes to 'slide down' 
the downhill part of the potential is 7(1 — a) 2 L 2 /Vo. 
Therefore, the discretization time has been chosen much 
smaller than these times, and also much smaller than 
the period of the rocking, 27r/f2, to avoid aliasing. In 
addition, the Langevin equations have been numerically 
solved up to times much greater than these characteristic 
times, and large enough to ensure a value of the com- 
puted flux close to its asymptotic value independently of 
the specific realization of the stochastic process, the ini- 
tial conditions, or the initial transitories. We have per- 
formed the number of realizations required to have small 
statistical errors in the averages computed. 



III. ONE-PARTICLE RATCHET 

For the one-particle ratchet (N = 1), the maximiza- 
tion of the instant velocity control policy, Eq. Q, only 
depends of the position x(t) of the particle. Hence, we 
can define an effective force F c g(x) = a(x)F(x) that al- 
lows us to rewrite the Langevin equation {1} as 



7±(f) = F eS (x(t)) + Acos(Q.t) + f (i). 



(6) 
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The effective force F e g derives from an effective poten- 
tial V c g(x) that is no longer periodic, but tilted down- 
hill. This V fi can be recast as a periodic potential 
Veff T ( x ) °f ne ight a ^o and asymmetry a, plus a linear 
term Vqx/L accounting for the bias, where Vb is the 
height of the ratchet potential, a its asymmetry pa- 
rameter, and L its period. Therefore, we can write 
V e g(x) = V^ T (x) — Vqx/L, as we illustrate in Fig. Q] for 
potential ([5]). In view of these considerations, the feed- 
back rocking ratchet can be reinterpreted as an open-loop 
rocking ratchet with a biased asymmetric potential. Thus 
Eq. © stands for the celebrated SQUID ratchet [ijj], 

im = ~V^(x(t)) + V /L + Aco S (nt) + £(t). (7) 

This equation of motion describes the dynamics of a 
tilted rocking ratchet, i.e., of a periodically driven sin- 
gle Brownian particle in a tilted washboard potential, 
and it has been extensively studied analytically and nu- 
merically [l2l. [l6l \Yl\ (even when inertial terms are also 
present [18|). For instance, for the adiabatic regime, i.e., 
the regime of slow driving [l3j , the flux can be approxi- 
mated by 




Jo 



where T = 2ir/Vt is the period of the driving force and 
(x)g(*) is the- asymptotic flux that would be obtained if 
the driving force were fixed at the instant t to its value 
G(t) = Acos(fit). This flux can be obtained by solv- 
ing a Fokker-Plank equation for the resulting constant 
external force || . Other analytical results have been re- 
ported for the high-frequency regim e |19| . or the deter- 
ministic (zero-temperature) regime [201 ] . Thanks to the 
equivalence found between the one-particle rocked feed- 
back ratchet [Eqs. {T} and @] and the SQUID ratchet 
[Eq. |(7J)], all the effects found for the rocked feedback 
ratchet have their counterparts in the extensively studied 
SQUID ratchet. However, it is important to emphasize 
that the tilt appears in our results not as an a priori bias, 
but as part of an effective description of the effects of the 
feedback. 

Here, we discuss the results for the different regimes 
obtained by performing numerical simulations of the 
Langevin equation. We shall first discuss the case of zero 
temperature and later the case of nonzero temperature. 

In the deterministic case of zero temperature there is 
no diffusion and only the rocking force can help the par- 
ticle to cross the flat regions of the effective potential V c s 
(see Fig. [T]). This makes the flux strictly zero for small 
amplitudes such that the particle cannot overcome the 
flat part of the effective tilted potential; see Fig. O For 
higher amplitudes the flux exhibits remarkable character- 
istic effects. Our simulations show that the deterministic 
flux is quantized and it presents a step-like structure, 
a well-known effect for open- loop rocking ratchets [l6l |. 
This structure is specially clear for the frequency Q = 50 
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FIG. 2: (Color online) One-particle case. Flux for the de- 
terministic (zero temperature) rocked feedback ratchet as a 
function of the amplitude A of the rocking and frequencies 
fl = 1,50,200 for the ratchet potential of Fig. [T] Inset: 
Zoom of the flux for Q — 50 and A belonging to the interval 
[48.0, 50.2]. Units: Vb = 5, L = 1, 7 = 1. 

in Fig. [2 The flux quantization is owing to the synchro- 
nization with the phase of the periodic driving (see [l6j 
for details). Its step- like structure presents a self-similar 
structure with steps at rational values of the flux, which 
can be seen performing successive zoom-in views (see in- 
set in Fig. 2). This structure is known as Devil's stair- 
case miii. 

Let us now discuss the case of nonzero temperature. 
A finite thermal noise leads to particle diffusion, which 
provides another mechanism to overcome the flat regions. 
This diffusion makes that the quantized step-like struc- 
ture for the flux is smeared and finally wiped out. On 
the other hand, a surprising effect is found for this case, 
namely a flux increase when the feedback policy and the 
rocking forcing arc both present. Unexpectedly, the re- 
sulting flux is greater than the sum of the flux values due 
to each separated effect, as we show in Fig. [3J In fact, 
the synchronization of the driving force with the feed- 
back mechanism gives positive large fluxes even for the 
case of negative fluxes for the pure rocking, i.e., with the 
ratchet potential always on [compare, for instance, curves 
for A = 40 in panels (a) and (b) of Fig. [31 or curves for 
fl = 100 in panels of Fig. [4]. Therefore, adding an exter- 
nal fluctuating force to the maximization of the instant 
velocity feedback protocol allows us to improve the per- 
formance of the system in a nontrivial way which to our 
knowledge has not been previously reported. This fact 
is not only relevant from a theoretical point of view, but 
also for experimental ratchet devices designed to maxi- 
mize the flux 0,0. 

Further insight on the behavior observed in panel (a) 
of Figs. [3J and 2] can be obtained studying the fast driv- 
ing regime. In this regime, it is useful to introduce a slow 
variable y(t) such that the position x(t) can be written 
as x{t) = y(t) + ip(t), where = rsin(f2t) is the fast 
contribution due to the fast driving, and r := A/("/Q). 
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When the driving is fast enough, a large number of oscil- 
lations in ip{t) take place before a significant change in 
y(t) occurs; thus we can proceed to the adiabatic elimi- 
nation of the fast variable tp(t) by averaging it over time. 
This procedure leads to an effective equation for the slow 
variable 



(a) 



7tf(*)=-FBff(!/(t)) +£(*), 
where F e s(y) = -V^ s (y), with 

1 



(9) 



V cS (y(t)) := - J V cS (y(t) + da. (10) 

This effective potential allows us to give a closed-form 
expression for the flux 



Lk B T 



1 _ e (VaniL)-V o!! (0))/k B T 



7 dx f* +L dy e(v-*(v)-vM*))/ 



k B T 



(11) 



Note that the potential V c s(y(t)) only depends on the 
characteristics of the driving force through the quotient 
r = A/("yQ,), and hence the same is true for the flux ob- 
tained within this fast driving regime. This approach is 
known as the vibrational mechanics scheme [21] . It has 
been successfully applied to the characterization of the 
so called vibrational resonance in bistable systems, both 
in the absence and presence of noise, as well as to the 
study of harmful effects (supression of the firing activ- 
ity) of strong, high-frequency fields on the response of 
excitable systems [22j]. In the context of ratchets, it has 
been used in the study of the effects of high frequency 
modulation on the output of Brownian particles moving 
in periodic one-dimensional substrates under the action 
of low- frequency input signals [23[ . The results obtained 
with this vibrational mechanics procedure are valid when 
the rocking force has frequencies much larger than the 
rest of characteristic frequencies of the system [2l[ . The 
average in Eq. ([TO)) makes the original potential barri- 
ers appear effectively lowered and flattened, eventually 
dissappearing as the ratio r increases. In particular, in 
our system, the periodic part of the one-particle effective 
potential V e g(x) = V^ T {x) — Vqx/L becomes smoother 
and smoother as an effect of the averaging process as r 
growths, and for large r only the linear term survives 
in the effective potential, giving a flux value of Vof(L'y). 
Panel (a) of Figs. [3] and 2] and Fig. [5] show how this value 
is reached for large amplitudes. Indeed, this is the largest 
value of the flux that has been obtained in a ratchet de- 
vice without an a priori bias; see Fig. [SJ 

The previous analysis also provides predictions on the 
dependency of the flux with the amplitude and frequency 
of the driving force. Within the vibrational regime, if the 
frequency is increased for a fixed amplitude, i.e., r is de- 
creased, then the flux will decrease until the value of the 
pure feedback ratchet [see panel (a) of Figs. [3] and [5]. 
Note that the values of the flux corresponding to low 
frequencies can be explained with the adiabatic descrip- 
tion in Eq. ©. On the other hand, if the amplitude is 
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FIG. 3: (Color online) One-particle case. Panel (a): Flux 
(x) versus frequency £1 for the rocked feedback ratchet. The 
horizontal solid line stands for the pure feedback ratchet — 
without rocking, i.e., A — 0. — Panel (b): Flux (x) versus 
frequency f2 for the pure rocking ratchet — without feedback 
flashing, i.e., a(t) — 1. — We have used the ratchet potential 
of Fig. [1] The one-sigma error bars are much smaller than 
the symbol size. Units: fcsT = 1, L = 1, 7 = 1. 



increased for a fixed frequency, i.e., r is increased, then 
the flux will increase from the value of the pure feedback 
ratchet up to the maximum value Vo/(L~/) [see panel (a) 
of Figs. Hand [5]. 

The effective potential V c g has allowed us to de- 
scribe the dynamics in the vibrational regime. We have 
compared the results obtained directly from Eqs. (fT0|) 
and (fTTj) with numerical simulations of the Langevin 
equation ([T]), with a good agreement for the fast driv- 
ing regime (Fig. [5]). This stress again the significance 
of the vibrational approach that has been revealed as 
a useful approach for both qualitative and quantitative 
predictions. On the other hand, we have found that 
fluxes greater than Vo/(Z>y) can be attained outside the 
vibrational regime in the quasideterministic regime, i.e., 
for large values of the potential height and the driving 
force amplitude. [For example, ratchet potential (JSJ) with 
Vq — 40, and a rocking force with A = 160 and f2 = 290 
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FIG. 4: (Color online) One-particle case. Panel (a): Flux 
(x) versus amplitude A for the rocked feedback ratchet. The 
pure feedback ratchet — without rocking — corresponds to 
the point A = 0. Panel (b): Flux (x) versus amplitude A 
for the pure rocking ratchet — without feedback flashing, i.e., 
a — 1. — We have used the ratchet potential of Fig. [1] The 
one-sigma error bars are much smaller than the symbol size. 
Units: k B T = 1, L = 1, 7 = 1. 

gives (i) ~ 43 in units fc^T = 1, L = 1, and 7 = 1]. 
This result is in accordance with the results found in 
Ref. [l7| for a tilted rocking ratchet in the quasideter- 
ministic regime. 

IV. COLLECTIVE RATCHET 

The dynamics of the collective ratchet compounded of 
more than one particle differs significatively from that of 
the one-particle ratchet discussed before. For collective 
closed-loop ratchets the feedback effectively couples the 
particles with each other and no simplifying description 
in terms of an effective potential has been found. 

The behavior of the deterministic (zero temperature) 
collective ratchet is similar to that of the one-particle 
ratchet, including the quantization of the flux and the 
step-like structure commented in Sec. IIIII We shall now 
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FIG. 5: (Color online) Flux (x) versus r — -^j for the one- 
particle rocked feedback ratchet. For increasing frequencies 
the values of the flux tend to the curve of the vibrational 
approximation. We have used ratchet potential (JSJ) with Vb = 
5. The one-sigma error bars are much smaller than the symbol 
size. Units: ksT = 1, L = 1, 7 = 1. 



focus in the nonzero temperature case for few- and many- 
particle collective ratchets where important differences 
emerge. 

In the few-particle case the maximum averaged center- 
of-mass flux is achieved for finite amplitudes and frequen- 
cies of the rocking force. Contrary to the one-particle 
case, the flux diminishes as the amplitude increases over 
its optimal value. On the other hand, we point out that 
for collective ratchets the maximum flux diminishes with 
the number of particles N. For a critical number of parti- 
cles the dependence of the flux with N practically disap- 
pears, indicating the transition to the many-particle case; 
see Fig. [S] The value for this iV-independent maximum 
flux that is obtained in the many-particle case coincides 
with the maximum flux obtained in the corresponding 
rocked flashing ratchet (open loop). This coincidence 
is analogous to the coincidence between the maximum 
flux for the threshold protocol in the many-particle case 
and the maximum flux obtained from the correspond- 
ing flashing ratchet [llj (see also Fig. [6]). Both of these 
coincidences can be interpreted as a consequence of the 
fact that these feedback protocols only use one bit of 
information about the system. This fact together with 
the increase of degrees of freedom of the system as N 
increases, makes that the relative strength of these feed- 
back protocols is weakened as the number of particles 
N increases, and that for systems with a large number 
of particles those feedback protocols cannot significantly 
beat their open-loop counterparts. In the following we 
discuss the interesting cooperative effects appearing in 
the many-particle case. 

In the many-particle case, the force per particle due 
to the ratchet potential, f(t) (defined in Sec. [TTJ) , has a 
quasideterministic evolution, as fluctuations in f(t) are 
subdominant. The analysis of f(t) has revealed to be 
very helpful to understand the dynamics. For the pure 
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FIG. 6: (Color online) Maximum center-of-mass flux versus 
number of particles N for the optimal rocked feedback ratchet 
(red +), the optimal threshold protocol (blue o), the instant 
maximization protocol — pure feedback — (green x), the op- 
timal rocked flashing ratchet (dotted line), the optimal peri- 
odic protocol (dashed line), and the optimal rocking ratchet 
— pure rocking — (solid line). We have used ratchet poten- 
tial (0 with Vo = 5. Where not shown, the error bars are 
smaller than the symbol size. Units: ksT = 1, L = 1, 7 = 1. 



feedback ratchet (without rocking) with many particles 
the system dynamics gets trapped with the potential 'on' 
or 'off' because the force fluctuations responsible of the 
switchings are negligible [4.]. Consequently, the system 
dynamics is near equilibrium most of the time and the 
net force is nearly zero. This implies an average asymp- 
totic center-of-mass velocity (xcm) tending to zero as N 
increases However, the introduction of the driving 
force allows the system to avoid this trapping and can 
result in an increase of the flux. 

Let us first discuss the cases of frequencies SI of lower 
or similar order to 2ir/Tf, with Tf the quasiperiod of 
f(t) for the pure feedback ratchet The maximum 
values of the flux in the many-particle case are obtained 
in this low-frequency regime. When the driving force is 
added, a complex synchronization appears between the 
quasideterministic dynamics of f(t) and the driving force 
Acos(Sli). We show in Fig. [7] [panel (a)] a typical time 
evolution of the forces for this case. The value of the 
flux depends on the details of this synchronization and 
it shows local maxima and minima when the system's 
parameters are tuned. See panel (b) of Fig. [71 where 
this complex behavior is shown by computing the flux 
for a two dimensional grid of 22 x 20 points in the A — ft 
plane. For the ratchet of this Figure the maximum flux 
(xcm) — 2.1 is achieved for a driving force of ampli- 
tude A ~ 20 and frequency SI ~ 55, expressed in units 
JcbT = 1, L = 1, 7 = 1. We want to call the attention 
to the fact that this frequency coincides with the charac- 
teristic frequency of the optimal threshold protocol [Tpj 
and of the optimal flashing ratchet SI = 27r/0.11 ~ 57. 
The maximum flux for the many-particle rocked feedback 
ratchet is reached when the rocking force has this charac- 
teristic frequency and pushes forward during the off pe- 
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FIG. 7: (Color online) Low and medium frequency rocking in 
the many-particle case (N = 10 4 ). Panel (a): Time evolution 
of the force / [Eq. <[3j ] for the rocked feedback ratchet (thick 
red line) and for the pure rocking ratchet (thin blue line), both 
with a rocking force (dashed line) of amplitude A = 10 and 
frequency SI = 20. Panel (b): Flux (xcm) versus amplitude 
and frequency for the rocked feedback ratchet. We have used 
ratchet potential ([5]) with Vo = 5. Units: fcsT = 1, L = 1, 
7=1. 



riod and backward during the on period. This makes that 
the ratchet potential only hinders the backward pushes 
of the rocking force. Contrary to the one-particle case, 
the flux diminishes as the amplitude increases over its op- 
timal value; compare Fig. [H panel (a), and Fig. [Jj panel 
(b). 

On the other hand, in the regime of large frequencies 
(SI 3> 2n/Tf) the pattern of fit) resembles the pattern 
for the pure feedback ratchet [J], but modulated by the 
high-frequency signal (Fig. [8|). For moderate values of 
the amplitude of the rocking, the system behaves more 
or less as if the fluctuations were increased. Therefore, 
an enlargement of the flux is possible for appropriate am- 
plitude of the driving force that succeeds in preventing 
the trapping similarly to the so-called threshold proto- 
col [10(. We show in Fig.[5]this resonantlike effect for this 
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averaged evolution equations: 
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FIG. 8: (Color online) Many-particle case (N = 10 5 ) high- 
frequency rocking (O = 1000) for ratchet potential ((5| with 
Vo — 5. Panel (a): Evolution of the force / [Eq. Q] for the 
feedback ratchet rocked with a high-frequency rocking force of 
amplitude A = 1 compared with the average force T,f=iFi/N 
with Fi given by Eq. (|13[) . We illustrate in the inset the mod- 
ulation of f(t) (thick red line) due to the high-frequency rock- 
ing (thin green line). Panel (b): Flux (xcm) versus A/(jQ) 
for two high-frequency rockings compared with the prediction 
of the vibrational approximation [Eqs. (|12p . (|13p ]. One-sigma 
error bars are shown. Units: fesT = 1, L = 1, 7 = 1. 



regime. We note that for small amplitudes A the system 
is not able to avoid trapping, while for too large ampli- 
tudes the characteristic quasideterministic f(t) pattern 
is erased and the flux goes to zero. 

As in the one-particle ratchet, a vibrational regime 
appears when the displacements induced by the driving 
force are faster than the effects of the other terms. This 
is shown in panel (b) of Fig. [5] In this panel, the de- 
pendence of the flux on the ratio r = A/(^Vt) for two 
different high-frequency rocking forces is compared with 
the flux obtained by assuming an effective dynamics de- 
fined as follows. As for the case N = 1, we introduce the 
slow variables yi(t) = Xi(t) — ip(t), with tp(t) — rsin(ilt) 
the displacements induced by the fast driving. Numeri- 
cal simulations confirm that the dynamics in this regime 
is governed by the slow variables verifying the following 



jy i {t) = F i (y 1 (t),...,y N (t))+Z i {t) ) 



(12) 



where 



Fi(yi,. 



T 



ds a(y 1 + ip(s), ...,y N + ip{s)) 

xFi(yi + ip(s), ...,y N + i>(s)), 

(13) 



with a given by Eq. (J4j) . This implies as before that 
within this regime the flux only depends on the charac- 
teristics of the driving force through the quotient r = 
A/(~fVt). We have numerically checked this for the few 
and the many-particle cases with high-frequency driv- 
ing forces, finding a better agreement in the few-particle 
case. However, we have also found a good agreement in 
the many-particle case for small values of the rocking am- 
plitude [see Fig.[8](b)] when £i=i F./N is a good average 
description of the force f(t) [see Fig. [5] (a)]. In addition 
to computations with ratchet potential ([5]), we have also 
performed computations with other potentials and found 
analogous results. 



V. CONCLUDING REMARKS 

In this paper we have studied the effects of rocking a 
feedback ratchet. The interplay between the rocking and 
the feedback policy gives an intriguing rich dynamics that 
we have analyzed and discussed. 

For the one-particle rocked feedback ratchet we have 
found an effective description in terms of a tilted rock- 
ing ratchet. Our simulations for rocked feedback ratch- 
ets show a relevant effect, namely, the magnification of 
the flux with respect to both the pure rocking and the 
pure feedback. That is, the rocked feedback ratchet is 
able to give fluxes even larger than the sum of the two 
fluxes separately. At this point, we remark that one of 
the main advantages of feedback ratchets over their open- 
loop counterparts is their ability to enlarge the particle 
flux, as it has been proved theoretically [4j and experi- 
mentally Q. In that sense, the introduction of the fluctu- 
ating force in feedback ratchets provides a way to further 
enhance the flux performance. In fact, the one-particle 
rocked feedback ratchet studied here gives the maximum 
flux that has been achieved in a ratchet device without 
an a priori bias (see Fig. [5]). This improvement in the 
flux performance is relevant for nanotechnological appli- 
cations of the ratchets. In addition, the observed depen- 
dence of the flux on the frequency and amplitude of the 
rocking signal has been explained for the whole range of 
parameters of interest. 

The rocking term also helps to enlarge the flux in the 
few- and many-particle case, as we have shown in Fig. [5] 
In this respect we highlight that the increase of the num- 
ber of particles effectively decreases the strength of the 
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feedback in the control, and in the limit of infinite number 
of particles this closed-loop protocol cannot give fluxes 
greater than its corresponding open-loop protocols. We 
have numerically shown the dependence of the flux with 
the amplitude and frequency of the driving force for these 
rocked feedback collective ratchets. The details of this 
dependence follow from the synchronization between the 
driving force and the feedback. In addition, we have 
found a new resonantlike effect when the amplitude of 
the rocking is tuned in the regime of high-frequency sig- 
nals. This later effect can be viewed as similar to an 
effective enlargement of the fluctuations in the net force, 
which prevents the trapping of the dynamics near equi- 
librium and results in an increase in the flux. 

To sum up, we have proposed and discussed a new 
closed-loop ratchet that is able to perform better than 
other known ratchets as a consequence of the nontrivial 
interplay of the feedback scheme and the rocking force. 
We have found an effective potential for the one-particle 



ratchet that explains the effective bias of the system and 
the a priori unexpected high values of the flux; it also 
has allow us to provide a closed expression for the flux 
in the vibrational regime. We have also analyzed the 
rich dynamics for the collective ratchet describing the 
vibrational regime and new resonantlike effects. 
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